Strongly-Correlated Photons Generated by Coupling a Three- or Four-level System to a Waveguide 
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We study the generation of strongly-correlated photons by coupling an atom to photonic quantum fields 
in a one-dimensional waveguide. Specifically, we consider a three-level or four-level system for the atom. 
Photon-photon bound-states emerge as a manifestation of the strong photon-photon correlation mediated by the 
atom. Effective repulsive or attractive interaction between photons can be produced, causing either suppressed 
multiphoton transmision (photon blockade) or enhanced multiphoton transmision (photon-induced tunneling). 
As a result, nonclassical light sources can be generated on demand by sending coherent states into the proposed 
system. We calcuate the second-order correlation function of the transmitted field and observe bunching and 
anti-bunching caused by the bound-states. Furthermore, we demonstrate that the proposed system can produce 
photon pairs with a high degree of spectral entanglement, which have a large capacity for carrying information 
and are important for large-alphabet quantum communication. 

PACS numbers: 42.50.Ct,42.50.Gy,42.79.Gn 



I. INTRODUCTION 

Strong coupling between light and atoms has been demon- 
strated both in classical cavity quantum electrodynamics 
(QED) systems and in more recent circuit-QED experi- 
ments ISHSl- This enables the generation of strong nonlinear 
photon-photon interactions at the single-photon level, which 
is of great interest for the observation of quantum nonlin- 
ear optical phenomena ll9l- [l2ll . the control of light quanta in 
quantum information protocols such as quantum networking 
lll3i [Till , as well as the study of str ongly c orrelated quan- 
tum many -body systems using light 115142311 . For example, 
both electromagnetically induced transparency (EIT) ifTlll and 
photon blockade 113, l25ll have been observed in recent ex- 
periments with trapped atoms in an optical cavity jH^-jl^ 
and with superconducting qubits in a microwave resonator 
ll29l [30I1 . Coherent transfer of quantum states between light 
and stationary qubits has been demonstrated in both cavity- 
QED jfsTl and circuit-QED if32l[33ll systems. In a very recent 
experiment, coherent transfer of photons between three res- 
onators has been realized in a superconducting circuit 1134)]. 

Recently, an alternative waveguide-based QED system 11351 - 
I47I1 has emerged as a promising candidate for achieving strong 
coupling between photons and atoms, motived by tremendous 
experimental progress ^ [13, l48l - l52ll . The experimental 
systems include a metallic nanowire coupled to a quantum 
dot lfT2ll . cold atoms trapped inside a hollow fiber ll48ll . a di- 
amond nanowire coupled to a quantum dot |49fl, a ID super- 
conducting (SC) transmission line coupled to a qubit |,8t, i29[l , 
and a GaAs photonic nanowire with embedded InAs quantum 
dots ifsollsill . In particular, it has been experimentally demon- 
strated that more than 90% of the spontaneously emitted light 
has been guided into the desired waveguide mode ifsill . deep 
into the strong-coupling 153(1 re g ime. Theoretically, single- 
photon switches 1381 |4U l44l 15411 have been proposed based 
on a waveguide QED scheme. An interesting photon-atom 
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oretically that EIT 11431 14511 and photon blockade 114511 emerge 
in a ID waveguide system. 

In this work, we consider using a waveguide-QED system 
to generate strongly-correlated photons through coupling to a 
three-level or four-level system (3LS or 4LS). Such strongly- 
correlated photons can be used to study many-body physics 
ifisll as well as to implement large-alphabet quantum com- 
munication protocols ll56[ IstIi . Specifically, to probe the 
strong photon-photon correlation mediated by the 3LS or 4LS, 
we study photonic transport, second-order correlation, and 
spectral entanglement of the correlat ed p hoton states. Fol- 
lowing Refs.|42, SH im and Issi and l59l we explicitly con- 
struct the scattering eigenstates by imposing an open bound- 
ary condition and setting the incident state to be a free plane 
wave. In the multiphoton solutions, photon-photon bound- 
states emerge, which have significant impact on the transport, 
spectral entanglement, and second-order correlation func- 
tion. While single-photon transport exhibits EIT, multiphoton 
transport shows photon-induced tunneling and photon block- 
ade. A highly entangled photon pair in frequency is obtained 
by scattering a two-photon state off the 4LS. Finally, we study 
the scattering of a coherent state wavepacket, whose number 
statistics become non-Poissonian. Strong bunching and anti- 
bunching appear in the second-order correlation function. 

This paper is organized as follows. In Sec. II, we intro- 
duce the model Hamiltonian, identify relevant experimental 
systems, and solve for the scattering eigenstates for one-, 
two- and three-photon states. With the scattering eigenstates, 
the asymptotic output states from scattering Fock states off 
the 3LS or 4LS are obtained in Sec. III. In Sec. IV, we study 
the photonic transport of Fock states and analyze the effect 
caused by the photon-photon bound-states. In Sec. V, we cal- 
culate the spectral entanglement for the two-photon case and 
demonstrate that highly entangled photon pairs are obtained. 
In Sec. VI, the signatures of photon correlation are revealed 
in the number statistics and second-order coiTelation function 
after scattering a coherent state wavepacket. Finally, we con- 
clude in Sec. VII. Some results related to photon blockade in 
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the 4LS were reported previoiusly in Ref.l45[ 



n. SYSTEM, HAMILTONIAN, AND SCATTERING 
EIGENSTATES 

We consider the scattering problem of photons in a one- 
dimensional waveguide side-coupled to a single atom, as 
shown in Figure[T] By "atom " we mean a local emitter with 
discrete levels, which could be formed from natural atoms, 
quantum dots, trapped ions, or superconducting qubits. 

Here, two types of local emitter are considered: a driven 
A-type 3LS and an A-type 4LS. The single-photon dynam- 
ics for the 3LS was previously studied in Ref.|4l] and a two- 
photon solution was found in Ref.|43|in the limit of weak con- 
trol field. Here, without assuming a weak control field, we 
solve the scattering problem for both the 3LS and 4LS in the 
general case. We mainly focus on the photon-photon correla- 
tion induced by the atom: physically, the interesting physics 
originates from the interplay of quantum intereference in the 
ID waveguide and interaction effects induced by the atom. 
Such interaction can be understood by treating the atom as 
a bosonic site and the ground and excited states as zero and 
one boson states, respectively. Unphysical multiple occupa- 
tion is removed by adding an infinitely large repulsive on- 
site interaction term [40], which is the underlying mechanism 
responsible for the formation of photon-photon bound states 
ill, m, m m m. The proposed system could be realized 
either in optical systems 1 12l pft dU l6ft 16111 or in microwave 
superconducting (SC) circuits Il8ll29l l62l - l64ll . For the optical 
systems, the driven 3LS and 4LS have been studied in both the 
trapped ion ll65ll and cavity systems ll27l l28ll66ll67ll . For the 
microwave SC systems, the 3LS and 4LS have already been 
reahzed using SC qubits IHIllllllillTl. 

We start with the Hamiltonian in the rotating wave approxi- 
mation, describing a continuurn photonic field in a ID waveg- 
uide coupled to a single atom ll36ll4ll - l43ll45ll 



H - H^o + //atom + He, 



dx{—i)hc 



+ a f a 

ag(x)—ai{(x) - fl[(x)— fli(x) 



(1) 



where ^(x) is the creation operator for a right- or left-going 
photon at position x and c is the group velocity of photons. 
For the driven A-type 3LS, 

^^om - Yj K'j-%J>(J\ + ^(|2X3| +h.c.), 

7=2,3 

H™ = /..«V«,)|[4,,).4,.),|l><2|.h.c.j. (2) 
For the A-type 4LS, 

<o!n = i:%-T)l^->0-|-?(|2><3Kh.c.), (3) 

7=2 

h['^^ = J dxnV6ix)[[al{x) + alixmi){2\ + 13)(4|) -Hh.c. 




(b) 





FIG. 1 . (color online) Sketch of the atom- waveguide system: (a) a A- 
type three-level system, (b) an A'-type four-level system, (c) photons 
(yellow) in a ID waveguide coupled to an atom (blue), which can 
be either the 3LS in (a) or the 4LS in (b). The transitions |1> «-> |2> 
and |3> «-> |4> are coupled to the waveguide modes with strength V. 
The transition |2> «-» |3) is driven by a semiclassical control field with 
Rabi frequency CI and detuning A. Here, oJc is the frequency of the 
control field. 



Here, the energy reference is the energy of the ground state 
and 62 = W21, £3 = £2- A, and £4 - £t, + W43, where W21 
and W43 are the |1) |2), and |3) |4) transition frequencies, 
respectively. In the spirit of the quantum jump picture iFtHi . 
we include an imaginary term in the energy level to model 
the spontaneous emission of the excited states at rate to 
modes other than the waveguide continuum.The spontaneous 
emission rate to the ID waveguide continuum is given by F = 
2V^/c (from Fermi's golden rule). Notice that the use of the 
rotating wave approximation is justified by the fact that HF «; 
hcL>2i, which is the case in current experiments lf8l l29ll50l - l52ll . 

It is convenient to transform the right/left modes to 
even/odd modes: al{x) - a^j^(x) + a'^{-x)l V2 and fl,^(x) = 

fl^(x) -aj^(-x)/ V2. This decomposes the Hamiltonian into 
two decoupled modes. The even mode couples to the atom 
and the odd mode is free: H = He + //« with 



(4a) 
(4b) 



r t d 

He- dx{-i)ncal,{x)—ae{x) + Ha_tom + Hi^^, 
J dx 

r t d 

Ho — I dx{-i)Tical,{x) — flo(x). 
J dx 

The coupling Hamiltonian He is now 

«.™=/„,*V*,|»l„)|l,<2|.hx.j, (5a, 

//W = ^ dxKVS{,x)[al{,x){\l){2\ + \'i){^) + h.c], (5b) 

where V - V2V. Hereafter, we will concentrate on solv- 
ing for the scattering eigenstates in the even space. Be- 
cause [//, he + «atom] - [H, ho] = for the number operators 
he/o = Jdx a'^i^J(x)aeio{x) and the atomic excitation «atom, the 
total number of excitations in both the even and odd spaces are 



3 



separately conserved. Therefore, a general n-excitation state 
in the even space (n = + natom) is given by 



J dx" g'^"\x)al(xi)---al(x„) (6a) 
^ i=23 

J dx" g^"\x) al{xi)---al{xn) 
- r dx"-'Y,ff^i^)S\jal{xi)---al{Xn.,) 



10, 1>, 
(6b) 



i=23 



-J dx--^ff\x)S\,al 



{xi)---al(Xn-2) 



10, 1>, 



where |0, 1) is the zero-photon state with the atom in the 
ground state |1) and S"^. - \j}{i\. 

The scattering eigenstates are constructed by imposing 
the open boundary condition that S^^(x) is a free-bosonic 
plane wave in the incident region ll42l l45l Issll . That is, for 

Xi , ■ • ■ ,Xn < 0, 



1 X— 1 e'^"' 

g^"\x) = — /lUi(XQi)---hlcMQn')^ l^kix) = — =, 



(7) 



where Q - {Qi,- ■ ■ ,Q„) is a permutation of (1, ■ ■ ■ ,«). Solving 
the Schrodinger equation with this open boundary condition, 
we find the scattering eigenstates for the systems we consider 
here (for a detailed derivation for a two-level system, see the 
Appendix of Ref.l42l). Below, we present the one-, two-, and 
three-photon scattering eigenstates, which have the same form 

I 



for the 3LS and 4LS cases. In the even space, the one-photon 
scattering eigenstate with eigenenergy E — Tick is given by 

g^^\x) = gdx) = hkix) [ei-x)+-tk9{x)] , (8a) 
[ck-£2 + A + iTi/2][ck - 62 + (/r2 - /r)/2] - ^2/4 



[ck-€2 + A + /Ts /2] [ck - 62 + (/r2 + iT)/2] 



(8b) 



where 6{x) is the step function. The one-photon scattering 
eigenstate is exactly the same for both the 3LS and 4LS be- 
cause it takes at least two quanta to excite level |4): for single- 
photon processes, the 3LS and 4LS cases are equivalent. 

For two-photon scattering, we start with a free plane wave 
in the region xi,X2 < 0, and use the Schrodinger equation to 
find the wave function first in the region xi < < X2 and then 
for < jci , X2 ll42ll . We arrive at the following two-photon scat- 
tering eigenstate with eigenenergy E - Tic{k\ + ^2): 



g^^\xi,X2) = ^y)^gkMQi)gk2(^Q2) 

Q 



(9a) 



PQ 



where f = {P\,P2) and Q - (Q\,Q2) are permutations of (1,2), 
9(xq.j) - 6(xq. - xqj), and B^^^ is a two-photon bound state — 
Re[yi,2] > 0. Our solution applies for the general case of ar- 
bitrary strength of the control field. Taking the weak control 
field limit for the 3LS case, we checked that one recovers the 
two-photon solution found in Ref.l43l 

Following the same procedure, we obtain the three-photon 
scattering eigenstate with eigenenergy E = Hc{ki +k2 + kj): 



g'^\xi ,X2,X3)^^lY,8ki (XQi )gk2 (XQ2)8k3 (XQ, ) + ^ [gkp^ {XQi )Bf^^ {XQ^ , ) 8(xq^) + Bf^ 1^^^ (XQ, , XQ^ , XQ3 ) 9(XQ^ )] k 

Bf^ A.^ (XQ, , XQ^ , XQ^ ) = e'^'"'^ +('^''2 +'<''3 '^-"03 ] [Z)i e-^' ^'Qi "-^ei I + D2 £"^2^23 "^e, I 

+D3 e->'il'^e3-^e2l-T2Ue2--^e,l + e-T2l-^'e3--^'e2l->'il^e2-'^eil]6»(xe3,)6'(xg2,), (10) 



where B*^^' is a three-photon bound state, P = {Pi,P2,P3) and 
Q = (61,62, Q3) are permutations of (1,2,3). The coefficients 
Ci,2 and £>i,2,3,4 in the bound states depend on the system 
parameters and have diiJerent functional forms for the 3LS 
and 4LS. Expressions for 712, Ci_2, and Di 2,3.4 are given in 
Appendix A. Notice that the bound states here have more 
structure than in the two-level case ll36[ IstI I42II : for exam- 
ple, the two-photon bound state has two characteristic binding 
strengths instead of one. This is due to the internal atomic 
structure: for the 3LS or 4LS, the photonic field couples to 
the transitions from the ground state to both of the eigenstates 
in the dressed state picture of levels |2) and |3), giving rise to 
two binding strengths. Such bound states are a manifestation 
of the photon-photon correlation induced by having more than 



two photons interact with the same atom. For the 4LS case, 
this leads to strikingly different multiphoton transport behav- 
ior compared to the single-photon transport ll45ll . 

From the scattering eigenstates, we construct n-photon (« = 
1 to 3) scattering matrices (S matrices) using the Lippmann- 
Schwinger formalism ll37ll42ll72ll . The output states are then 
obtained by applying the S matrices on the incident states ll42ll . 

III. OUTPUT STATES OF FOCK STATE SCATTERING 

In this section, we present the output states from scattering 
one-, two-, and three-photon number states off of a 3LS or 
4LS. We assume that the incident state propagates to the right 
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and the atom is initially in the ground state. Specifically, we 
consider incident states in the form of a wavepacket for two 
reasons: (i) in practice, any state that contains a finite number 
of photons is a wavepacket; (ii) as we will show, sending in 
wavepackets with a finite width is crucial in order to observe 
the bound state effects in the measurements. The continuous- 
mode photon-wavepacket creation operator is given by iFtsIi 



■ Jdk Q-i 



(11) 



where a'j^iJJc) = (1/ V27r)/ dx e'^''^aj^^^(x) and the amplitude 

aik) satisfies the normalization condition ^ dk \a(k)\^ - 1. An 
incident right-going «-photon Fock state is defined as 



\na)R = 



(a' Y 



|0>. 



(12) 



With the n-photon S matrices S^'^\ we are able to find the 
asymptotic output state long after the scattering (f — > +oo) 
14211 . Specifically, the single-photon output state is given by 



(13a) 

(13b) 
(13c) 
(13d) 



|(/r(')) = J dka{k)\(f>^^\k)), 

\<P^'\k)}^tk\k}R + rk\k}L, 
\k}R/L = al^^(k)m, 
tk = (ik + l)l2, n = {tk-l)l2. 
The two-photon output state reads 

|(A*2)> = ^ dkidk2-^a(h)a(k2)\4>'^^\ki,k2)), (14a 
\(p^'^\ki,k2)) = J" dx\dx2\)^tki,k2{x\,X2)a\(xi)a\{x2) 



+rtk, ,k2 (xi , -X2)a]^{xi)a[{x2) 

1 t t 1 

+ 1^^, mJ-xi , -X2)fl^(xi )a^(x2)J|l 



(14b) 



where 



1 



tkt,k2 = tkitkihkM^y^kiixi) + ^Bflk^(xi,X2) + ki <-»^2 



1 



rtk.M = tk, rkJu. (xi)/z<;,(x2) + ^Bf\(xi,X2) + ki <^ k2 



4 ki,k2^ 

-B?\ V 

4 ki,k2^ 



rki,k2 = rkirk2hkMi)'^k2ix2) + ^B\^ j^^{xi,X2) + ki ^ k2, 

3(2) 

^\x\,A2) = t: ' - - 

j=l,2 



B'^X^xuX2)^e'^'^^'^^'^ 2 C,^^"''^l-^'2-^il0(x2i) + (xi ^X2). 



(15) 



In Eq. ( fT4l l. the output state has three components tk^ ,At , 
ftkiM (which is not a product), and rk^^k2^ corresponding to 
two-photon transmission, one-photon transmitted and one- 
photon reflected, and two-photon reflection, respectively. The 
first term in each of these functions is the plane-wave term. 
The second term is the bound-state term associated with the 
momentum-nonconserved (for individual photons) processes. 



(a) a = 0.01,Q = (b) a = 0.01, Q = 1.6 




(c) a = 0.2, a = 




1 




4-2024 

(d) a = 0.2, n = 1.6 


1 

— / 


1 



FIG. 2. Single-photon transmission T (solid), reflection R (dashed) 
and loss (dotted) as a function of incident photon detuning, for the 
values of cr (the wavepacket width) and £1 (the strength of the control 
field) shown. Here, the effective Purcell factor is P = 9. Note the 
sharp EIT window, particularly in the narrow wavepacket case. 



The three-photon output state takes a similar form and is 
shown in Appendix B. 

With the output states, we can study induced photon-photon 
correlation by applying various measurements on them. We 
present results for transport, spectral entanglement, number 
statistics, and second-order correlation in the following three 
sections. Throughout the paper, we choose incident Gaussian 
wavepackets with the spectral ampHtude 



a(cij) ■■ 



1 



(27ro-2)l/4 



exp 



4o-2 



(16) 



where cr is the width and ojq is the central frequency. We 
assume that level |3) is metastable (r3=0) and levels |2) and 
|4) have the same loss rate: r2 = In addition, we assume 
that the transitions |1) |2) and |3) «-» |4) are at the same fre- 
quency, W2I - W43, and the detuning of the control field is 
zero, A = 0. We set the loss rate as our reference frequency 
unit: r2 = r4 = I. The coupling strength to the waveguide 
is characterized by the effective Purcell factor P = r/r2 = F. 
Plasmonic waveguide swtems have been predicted to have 
a large Purcell factor OSll and a value of P = 1.5 has been 
demonstrated experimentally lfT2ll . Slot waveguides have been 
theoretically shown to have large values of P reaching 16. Re- 
cently, by carefully tailoring the ends of photonic nanowires, 
J. Claudon et al. achieved a value of f > 9 in the experiment 
ifsol Isill . Furthermore, 5.7 < f < 24 was demonstrated in a 
photonic crystal waveguide coupled to a quantum dot ll74ll . In 
superconductirig circuits with ID open superconducting trans- 
mission lines 13, l52l| . even larger values of P have been 
achieved, exceeding 15 ll52ll . 
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FIG. 3. (color online) Two-photon transmission and reflection probabilities for the 3LS and 4LS cases, (a)-(c) As a function of incident photon 
detuning 6at with P = 9 and cr = 0.2. (a) Probability that both photons are transmitted (and hence are right-going, P^l^)- (b) Probability that 

one photon is transmitted and one reflected (right-left, P^^l)- (c) Probability that both photons are reflected (both left-going, P^H)- (d)-(f) As 
a function of P with Sai = and cr = 0.2. (g)-(i) As a function of cr with P = 9 and 6ai = 0. The label PW refers to the contribution from the 
plane-wave term only, while BS refers to all the other contributions involving bound-state terms [Eq. l|19ll1. Here, we set Q = 1.6. The bound 
state effect enhances transparency in the 3LS case but blocks two-photon transmission past a 4LS. Note that a non-zero cr is crucial to observe 
these effects. 



IV. TRANSPORT OF FEW-PHOTON STATES 



B. Two-Photon 



A. Single-Photon 

With the output state in Eq. ( fT3l ), the transmission (T), and 
reflection (R) probabiHties for a single-photon are 

r = J dk\R{k\i(,'-^^}f^ J dkaHkm\ (17a) 



dk\L{k\l/,'-'^}f . 



I 



dka^(k)\rkf-. 



(17b) 



which are the same for both the 3LS and 4LS cases. Fig- 
ure|2] shows T, R, and the loss {I - T -R) as a function of 
the detuning 6co = <^o- <^2i at P -9. Clearly, EIT appears in 
Fig Elb), when the control field is on. As one increases the 
width of the wavepacket, as shown in Fig.|2jd), the BIT peak 
is suppressed as cr becomes comparable with the width of EIT 
window (~ Q^/F), see Eqs. ([8]) and ( fT3] ). In Fig.lSta), and (c), 
we set Q = 0, which means the control field is ofi" and the 3LS 
(4LS) becomes a reflective two-level system SlUlil]. No- 
tice that the width of the reflective peak in the Q.-Q case is 
~ F and hence is insensitive to the increase of cr from 0.01 to 
0.2. 



The two-photon transmission and reflection probabilities 
are given by 



pL^)^ I dkidk2^\RR{ki,k2\<P^^Y, 



RR 

0(2) 
RL 



dkidk2\RL{kuk2\>/^^^^}\\ 



Pfi^ j'dkidk2^\LL{kuk2\'/^'^^^}\\ 



(18a) 
(18b) 
(18c) 



where P 



.(2) p(2) 



(2) 

^ and P^^ are the probabilities to observe two 
transmitted photons, one transmitted and one reflected pho- 
tons, and two reflected photons, respectively. We separate the 
two-photon transmission and reflection probabilities into two 
parts: (P^^')pw is the contribution from indepedent single- 
particle transmission (denoted PW for "plane wave"), and 
(P*'*)bs is the contribution from both the bound-state term 
in Eq. ( fT4] i and the interference between the plane wave and 
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(2) 

bound-State terms. As an example, f is split as follows 



PrR = j dhdk2\Uh,k2) + B{ki,k2t 
(OpW = J dkidk2\t2{kuk2)\\ 

(PrXs^ j dhdk2[r^{kx,k2)B{kuk2) 

+t2{kuk2)B\kuk2) + \B{kuk2)\\ (19c) 



(19a) 
(19b) 



where 



hikiM) = a(ki)a(k2)tkitk2, 

B{kiM) = T- y ii-^ + 1-^ — ) 



(20) 



7=1,2 



; ^ dk a(k) aiky +k2-k) Cjik, ki+k2-k). 



Figure[3] shows the two-photon transmission and reflection 
probabilities for both the 3LS and 4LS cases, decomposed in 
this way. Because the PW term is from the single-particle 
solution, it is the same for both the 3LS and 4LS. However, 
(P<2))gs is quite different for the 3LS and 4LS. Figure[3l;a)-(c) 
shows F'-^' as a function of incident photon detuning. Close to 
resonance, in the 3LS case (P*^^)bs enhances the two-photon 



transmission P^^j^ while suppressing P^^l^. In contrast, in the 
4LS case (^'^^^)bs has exactly the opposite effect. This leads to 
enhanced multiphoton EITfor the 3LS i^il/ and photon block- 
ade for the 4LS M^- Such enhanced EIT and photon blockade 
are caused by the interference between the two multiphoton 
scattering pathways: passing by the atom as independent par- 
ticles or a composite particle in the form of bound states (for a 
detailed analysis, see the Supplementary Material of Ref.l45l). 
In Fig.[ld)-(f), we plot P^^) 

as a function the effective Pur- 
cell factor P for the on-resonance case. So, - 0. It is remark- 
able that, for f and f j^^, (^'*^^'')bs becomes comparable to 
(P*^^)pw in the strong coupling regime. An important impli- 
cation is that the bound-state effect can be observed in pho- 
tonic transport experiments, given recent rapid experimental 
advances fllHEBSl- 
Fig.[3g)-(i) shows 

as a function of the wavepacket 
width cr with P - 9 and the photons on resonance with the 
atom. There are several notable features. First, as cr ap- 
proaches zero, (P^'^^)-B5 shrinks to zero for both the 3LS 
and 4LS cases. This further highlights that sending in a 
wavepacket with a finite width is crucial to observe the bound 
state effect in photonic transport. Physically, this occurs be- 
cause, in the cr = limit under BIT conditions, the atom is 
fully transparent {T - 1) to the incoming photons and hence 
the atom-mediated photon-photon interaction is absent, in- 
hibiting any bound state effect. For the general case without 
EIT conditions, the above conclusion still holds: as cr —> 0, 
the bound state effect vanishes in multiphoton transport. This 
is because the bound-state term in Eq. ( l20b originates from 



.(2) 



3LS 

10'* E-(a) P = 9,(J = 0.2 



4LS 

(d) P = 9.(T = 0,2 




FIG. 4. Photon blockade and photon-induced tunneling in transmis- 
sion. Photon blockade strengths P21 (solid line) and P31 (dashed 
line) as a function of incident photon detuning 6oj, P and cr for (a)- 
(c) the 3LS case, and (d)-(f) the 4LS case. Here, = 1.6. The 3LS 
causes photon-induced tunneling while the 4LS causes photon block- 
ade. 



the coincident photons at the atomic site: as cr ^ 0, the 
wavepacket becomes infinitely long and the probability of co- 
incidence vanishes. Second, notice that while (P^^*)bs ap- 

(2) (2) 

proaches zero for /"^^^ as cr increases, its magnitude for f 

and P£l^ increases after an initial decrease. This is due to 
the enhanced interference between the plane-wave and bound- 
state terms [Eq. ^] for P^^l and pf/. 

The result for three photon scattering shows behavior sim- 
ilar to the two-photon case. To avoid duplication, we do not 
present it here. 



C. Photon Blockade and Photon-Induced Tiinneling 

To quantify the observed enhancement of FIT and photon 
blockade in Fig. [3] we define the strength of photon block- 
ade P21 for the two-photon case by the conditional probabil- 
ity for transmitting a second photon given that the first photon 
has already been transmitted, normalized by the single-photon 
transmission probability. Similarly, we can define P31 for the 
three-photon case. We thus have 



p(2) 
p - RR 



P31 = 



d(3) 
RRR 

r3 



(21) 



As shown in Fig.|4ta)-(c), for the 3LS case, the single-photon 
FIT is enhanced in two-photon and three-photon transmission 
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OJl-UJo LOi-LUq UJi-UjQ 



FIG. 5. Two-photon joint spectrum of the output states after scattering off a 3LS and 4LS in the case of a spectrally narrow incident wavepacket. 
Panels (a)-(c) show the uncorrelated spectra T2{cl)i,cl)2), RT((jj\,ciJ2), and R2{(jj 1,0^2)^ respectively. Panels (d)-(f) show the joint spectra of the 
transmitted field after scattering off a 3LS for two transmitted photons, for one transmitted and one reflected, and for two reflected photons, 
respectively. Panels (g)-(i) show the joint spectra of the transmitted field after scattering off a 4LS. Strong spectral entanglement is indicated 
in panel (i); this reflected field is essentially a pure two-photon bound state. System paramters: P = 9, fl = 1.6, (5a) = 0, cr = 0.01. 



by interaction with the 3LS. Pronounced photon-induced tun- 
neling iItsIi due to the strong correlations between transmit- 
ted photons occurs in this case: P2i,Pii > 1. In contrast, 
as shown in Fig.|4jd) and (e), scattering from a 4LS exhibits 
a different behavior within the BIT window, namely, photon 
blockade 1I45II : P2i,P3i < 1- For increasing coupling strength 
[Fig. lite)], P21 and P31 approach zero asymptotically when 
the incident photons are on resonance with the 4LS. In ad- 
dition, from Fig.|3c) and (f), we confirm that both photon 
blockade and photon-induced tunneling go away in the zero- 
width limit (cr — > 0). 



V. SPECTRAL ENTANGLEMENT OF PHOTON PAIRS 

It is clear that the two-photon bound state in Eq. ( fT4l i is 
entangled in the momentum (or equivalently frequency) de- 
gree of freedom. To probe this spectral aspect of the two- 
photon entanglement, we rewrite the two-photon output state 
[Eq. (fT4l il in frequency space as 

I'A'^^') = J dLi)idLi)2[fRR{Li)uCi)2)al,((iJi)al{ciJ2) (22) 

+fLL{(^i , W2)a[('^l)fl[(<W2)]|0), 

where //;/j(wi, ^2), /ffL(wi,W2), and /ll(wi,W2) are the two- 
photon amplitudes for a transmitted pair, a pair of one trans- 
mitted and one reflected, and a reflected pair, respectively. Ex- 



plicitly, they take the following form 



fRR{aJuU2) 


= f2(wi , W2) + B{loi , (02), 


(23a) 


fRdOJ\, 0)2) 


- 2[rt(a>i,LL>2) + B(aji,a>2)], 


(23b) 


fLL{0i\,(^2) 


— r2(wi,t<J2) + B{(jJi,(jJ2), 


(23c) 


?2(W1,W2) 


- ta,itoj20'{aJi)a(a>2), 


(23d) 


ft{cj\,0J2) 


= tajiraj2a((Dl)a((jJ2), 


(23e) 


r2(a)i,W2) 


= raj,ra,2a(<jJi)a((D2), 


(23f) 



where B(aJi,W2) is given in Eq. (|2Q]i. The first term in 
f{a>i,a>2) is the uncorrelated contribution, while the second 
term signals photon correlation. From Eq. (l23T l. we define the 
joint spectral function of the two-photon states to be iItgIi 

Fafi=RR.RL.LL(m,(^2) = l/a-/j(wi,W2)l^ • (24) 

For the purpose of comparison, we also define the uncorre- 
lated spectral function of the two-photon states, 

r2(wi,W2) = |f2(wi,W2)P, (25a) 

RTioJi,uJ2) = 4\Ft(a)uoJ2)\^, (25b) 

R2(cOl,CLl2) = \f2{ui,U2)\^ ■ (25c) 

Note we are using RT to denote the joint spectral function for 
one transmitted and one reflected (uncorrelated) pair. 

Figure|5]shows the two-photon uncoiTelated and joint spec- 
tra in the case of on-resonance photons (6a, - 0) and for a spec- 
trally narrow wavepacket (cr - 0.01). With the chosen param- 
eters, the FIT peak width is much larger than the wavepacket, 
~ O^/r = 0.28 » cr. Therefore, for the uncorrelated pair of 
transmitted photons [T2, Fig. 13 a)], there is only a sharp peak 
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FIG. 6. Nonclassical light source. Photon number statistics quantified by logio(/'„//'„ poi,sson)j where P,, and Pn,Poisson are the )2-photon 
probability in the transmitted field and in a coherent state with the same mean photon number, respectively. Panels (a)-(d) show the results 
of the transmitted field after scattering off the 3LS for n = 0, 1,2,3, respectively. Panels (e)-(h) show the results of the transmitted field after 
scattering off the 4LS for « = 0, 1,2, 3, respectively. The dashed line is a guide to the eye for equal probabilities, -Pn/^n.Poisson = 1- The speckle 
in the plots is numerical noise, coming from numerical evaluation of high-dimensional integrals in computing the transmission and reflection 
probabilities. System paramters: P = 9, Q.= 1.6, = 0, cr = 0.01, and 77 = 1 in the incident coherent state. Scattering off a 3LS enhances 
the multiphoton content of the pulse because of multi-photon EIT; in contrast, the photon blockade in the 4LS case suppresses essentially all 
multi-photon content, thus realizing a single-photon source. 



at LOi - a>2 - cl>o caused by the Gaussian spectrum of the in- 
cident photons. For the uncorrelated pair of one transmitted 
and one reflected photons (RT), there are two peaks resulting 
from the interplay of the spectrum of the incident photons and 
the rapid increase of the reflection probability away from the 
EIT peak (see Fig.|2|. Accordingly, there are four peaks for 
the case of two reflected photons, as shown in Fig.|5lc). 

Figure|3d)-(f) shows the joint spectra for the case of 3LS 
scattering. It is evident that the joint spectra of the pair of 
two transmitted photons {Frr), and the pair of one transmit- 
ted and one reflected photons (Frl), are dominated by the un- 
correlated transmission. The joint spectrum of the pair of two 
reflected photons [Fig. 13 f)] is slightly modified from the un- 
corrected spectrum along the diagonal line. This is caused by 
the coiTelated bound state term B{cji,a)2). For the 3LS case 
with the chosen parameters, the correlation term B{a)i,a>2) is 
of order 10~' and hence is too weak to affect Frr and Frl. 

In contrast, for the 4LS case [Fig.|5jg)-(i)], Fri and Fn 
are greatly modified by the correlation term, while Frr is still 
dominated by the uncorrelated transmission. In particular, as 
shown in Fig.|5](i), the joint spectrum of the reflected pair is 
dominated by B(a)i,a>2)- This pair is primarily made up of a 
pure two-photon bound state: the frequencies of the photon 
pair are correlated along the line coi +cl>2 - 2cl>o with uncer- 
tainty cr. A similar correlated photon pair was obtained in a 
waveguide-cavity system fl%. 

The two-photon bound state is a composite object of pho- 
tons with effective attractive interaction; it displays strong 
bunching behavior in photon-photon correlation measure- 
ments. Such a photon pair is highly entangled in frequency 
because measurement of the frequency of one photon unam- 
biguously determines that of the other. This strong spectral 



correlation provides more information per photon pair and 
could be used to implement large-alphabet quantum commu- 
nication IstIi . 



VI. COHERENT-STATE SCATTERING 

In this section, we study the scattering of a coherent state off 
a 3LS or 4LS. We probe the strong photon-photon correlation 
in the transmitted field by studying first the number statistics 
and then the second-order correlation function. 



A. Number Statistics 

We consider the case that the 3LS or 4LS is in its ground 
state initially and there is an incident continuous-mode coher- 
ent state of mean photon number « = 1, spectral width cr - 0.2, 
and central frequency on resonance with the atom, oiq - <jD21- 
In this case, the contribution from the four-photon state can 
be neglected (~ 1.6%). The photon-number statistics in the 
transmitted field is obtained by first applying the S matrices to 
the incident state and then measuring the transmitted field, as 
described in Ref.l42l 

We present the results for both the 3LS and 4LS cases in 
Fig.|6]by taking the ratio of the photon-number distribution in 
the transmitted field f „ (n = 0, 1 , 2, 3) to that of a coherent state 
^'«,Poisson having the same mean photon number as the trans- 
mitted field. From Fig.|6|a)-(d), it is clear that when the FIT 
condition is satisfied, the 3LS induces strong photon-photon 
interactions, which in turn reduce the one-photon probability 
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5 10 15 20 

r 

FIG. 7. Bunching and anti-bunching. Second-order correlation func- 
tion g*^^'(T) of tfie transmitted field for a weak incident coherent state 
(n <gi 1) of width cr = 0.2, resonant with the atom {6aj = 0). (a) Color 
map plot of logio[g'^^'(0)] as a function of the strength of the classical 
control field, CI, and the effective Purcell factor P. The dashed line 
marks the border between bunching (g*^\0) > 1) and anti-bunching 
(g*^'(0) < 1) behavior, (b) g'^* as a function of time delay t in four 
cases (using Q = 1.6). Inset: zoom at short time scales. 



and redistributes the weight to the two- and three-photon prob- 
abilities. This comes about because the bound state in the 3LS 
case enhances multiphoton EIT, as we have shown in Sec. IV 
B and C. 

In contrast, for the 4LS case shown in Fig.|6je)-(h), in most 
of the parameter space, we have enhanced single-photon prob- 
ability while suppressed multiphoton content: Pi > Pi,poisson 
and -P2(3) < P7.n\Pniss nii. This gives rise to a sub-Poissonian 
single-photon source f45], which comes about because, while 
EIT occurs in the single-photon transmission, multiphoton 
states experience photon blockade, as shown in Sec. IV B and 
C. Therefore, we demonstrate that the waveguide-atom sys- 
tem is capable of generating nonclassical li ght, which may 
find applications in quantum cryptography IITSUSlll or dis- 
tributed quantum networking |[l3l UM . 



B. Second-order Correlation 

To further probe the nonclassical character of the transmit- 
ted field, we calculate the second-order correlation function 
g^^^r), which is often measured experimentally. For a steady 
state, of the transmitted field is defined as 

(aL(x,t) aL(x,t + t) ap(x,f + t) aj/ixj)) 

g(2)(T) = Iim— ^ . (26) 

'^°°<fljj(x,f) aR{x,t)){a^j^{x,t + t) aR{x,t + T)) 

As shown in Appendix C, for our system, this definition is 
equivalent to following expression in the Schrodinger picture, 

<i/^|fl^(x) aR{x)\il/){ip\a\{x + ct) agix + ct)\i//} 

where Itfr) is the asymptotic output state. With a weak inci- 
dent coherent state (mean photon number n <«; 1), we con- 
sider only the contribution of the two-photon and one-photon 
states in the numerator and denominator in Eq. (IZTT i. respec- 
tively. Substitution of the single-photon and two-photon trans- 
mission wavefunctions from Eqs. ( fT3l l and ( fT4l i into Eq. (ITTb 
yields the explicit expression 

I / dkidk2 aih) a(k2) [fA, tk^ie-''^' + e-''^ + B(r)]p 

g (t) - — -— , 

I Jdkidki aih) aiki) tk,tk^(e-'''i^ + e-^'<2r)\2 

B{t) = n{Cie-^'" + Cie'^'^") . (28) 

In the numerator, the first term and the second term B{t) come 
from the plane wave and bound state pieces, respectively, in 
Eq.ldli. 

Figure|7la) shows §*^'(0), which is the same for the 3LS 
and 4LS cases. The presence of level |4) does not contribute 
to g^^\Q): it takes two quanta to excite |4), which then deex- 
ites in the form of cascade emission with zero probability to 
emit two photons at the same time. In Fig.|7la), there is rich 
bunching and anti-bunching behavior, caused by the two-body 
bound state. At r=0, the amplitude of the bound state term in 
Eq. (|28]l is Z?(0) = -Iri^rk^, where ri^^^ is the single-photon 
reflection coefficient. Hence, in the numerator of g^^'(O), the 
amplitudes of the plane-wave and bound-state terms are out 
of phase. When P = 0, the bound state term is zero and 
g^^\Q) - I. As, P increases, the strength of the bound state in- 
creases, causing ^^^'(0) to decrease until the bound state term 
cancels the plane wave term exactly, producing complete anti- 
bunching. Further increase of P leads to a rise of §(2)(0) and 
eventually photon bunching. 

By comparing Fig.|4lc)-(d) and Fig.|2ta), we find that pho- 
ton anti-bunching and photon blockade, and photon bunching 
and photon-induced tunneling do not have a one-to-one cor- 
respondence. For example, in the whole parameter regime of 
of Fig. Hid), photon blockade is present; while in Fig.|7] there 
is a large region of parameter space where photon bunching 
[g^^\Q) > 1] instead of photon anti-bunching [g^^\Q) < 1] is 
observed. This is because we are studying a state of contin- 
uous modes and performing instantaneous measurements at 
two space-time points (x, f) and (x, t + r). If one integrates over 
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the time t in the measurement iItsIi . as done in many experi- 
ments in which the detector integration time is much longer 
than the wavepacket duration, one finds a one-to-one corre- 
spondence between photon anti-bunching and photon block- 
ade, and photon bunching and photon-induced tunneling. 

The time dependence of g^^^r) is shown in Fig.|7lb). There 
are two characteristic time scales: tj - 1/Re[cyi] and T2 - 
1 /Re[cy2]- Within the short time scale, g*^* can display either 
bunching or anti-bunching for both the 3LS and 4LS cases, 
depending on the system parameters, as shown in the inset of 
Fig-H On the long time scale, for the 3LS case, g'^' shows 
bunching — g^^Hr) > 1 — corresponding to the enhanced mul- 
tiphoton transmission akeady apparent from both the photon- 
induced tunneling [Fig-Sc)] and the enhanced multiphoton 
content in the number statistics [Fig.|6l. For the 4LS case, 
anti-bunching [§^^'(t) < 1] dominates at long times, corre- 
sponding to the photon blockade observed in Fig.Ud) and 
the enhanced single-photon content in Fig.|6l Hence, for our 
pulsed output state, g^^\T = 0) displays rich physics due to 
the induced photon-photon correlation, but is not necessarily 
a good guide to the photon statistics. 

VII. CONCLUSIONS 

In summary, we present a waveguide-QED-based scheme 
to generate strongly-correlated photons, of interest for both 
many-body physics and quantum information science. Pho- 
ton bound-states appear in the scattering eigenstates as a man- 
ifestation of the photon-photon correlation. As a result, while 

I 



a single-photon experiences BIT in the proposed waveguide- 
atom system, multiphoton states can display either photon 
blockade or photon-induced tunneling, depending on the de- 
tailed structure of the "atom". From either the photon block- 
ade or photon-induced tunneling that occurs, nonclassical 
light sources can be generated by sending coherent states into 
the system. In the most interesting case, a 4LS removes the 
multiphoton content from the coherent state, leaving a pulse 
with only zero or single photon content. 

In addition, we find that the system can be used to produce 
highly entangled photon pair states in frequency space, po- 
tentially of use for large alphabet quantum communication. 
Finally, we show that rich bunching or anti-bunching behav- 
ior is present in the second-order correlation function as a 
signature of the strong photon-photon correlated mediated by 
the "atom". Given the recent rapid experimental advances in 
several realizations, the proposed waveguide-QED system is 
emerging as a promising route to cavity-free open quantum 
networks, which are crucial for both large-scale quantum com- 
putation and long-distance quantum communication. 
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APPENDLX A: EXPRESSIONS FOR 71,2, Ci 2 AND £'1,2,3,4 

In this Appendix, we give explicit expressions for the constants yi 2, Ci,2, and £>i,2.3,4 that appear in Eqs. (9) and (10) for 
both the 3LS and 4LS scattering eigenstates. yi,2 is the same for both cases and is given by 



(Ala) 
(Alb) 



FH-Fo+Fs /A \ rH-F2-Hr3 /A \ 

cyi = ^+iy- + e2 + riy cyi^ +^-iy--e2-riy 

For the A-type 3LS and A'-type 4LS cases, Ci,2 and 01,2,3,4 take the same form 

riK,N),, I^'''''\kuk2)-a{kuk2)^2 ^(,a,n),. -f^'''''\kuk2) + a{kuk2)Ai 

C '(^1,^2) = , C\' '{ki,k2)^ ^ , 

A1-A2 ^ A1-A2 

DrWk2,k,) . <^V^0-'^.3(^.)^. ^,,,,,^^^^^^^^ D<^.-)(.a2,^3) = "^^^'^'^''^7"^''^"' cf->(.2,^3), 

A1-A2 ' ^ A1-A2 ^ 

Dt-\kM.) - "^■^'''^'^^7"^'^^^^ cf'-'(^.^3), Dr\kM.) - ^^r^(^o-a.4(^i)^. ^,.,,>^^^^^^^^ 

A\-A2 ^ ^ A1-A2 

rH-F2-r3 /A \ r+F7-r3 /a \ 

^1 = 1 '-+^ + i[- + Ti], ^2 = 1 i-^+/(--77), (A2) 
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where the superscript A stands for the 3LS and for the 4LS. a's and yS's in the above equation read 



a(ki,k2) - — 



(A) . 



v{ki,k2) = 



In 

h, - 1 Ik, 
— + — 



-1 



Pi-2 Pki 

e4-£-(/r4-iT)/2 



(N) 



~tki-v{k\,k2) ~tk,-v{ki,k2) 



Pk. 



/T3 



e4-£'-(/r4 + /r)/2' 
where 7<. is given in Eq. ([Sj)) in the main text, a^, 024, j6i3 and ^624 are given by 

2{lk-l) 



Pk^\ck-£2 + ^+ — \\ck-e2+ 2 



Pki 

ir2 + /r 



Q2 

4 ' 



anik) = a24{k) ■■ 



V2^ 



f3i3(ki^^ = 
/324(kY^^ = 



1 



V2^ 
1 



V2^ 



4pk 
4pk 



-(fA-l)^l 



/3,,ikr>^^l^-\i,-^^mU 

V2^ 1 4pi L J 

-[ri:-jU2W] A2 



P24{kr- 



1 fm^ 



V2^ l4pi 



64 — /r4/2 — cA; + /r/2 + /cyi^ 
64 - ir4 /2 - cA - /r/2 + icy 1 ,2 



(A3a) 
(A3b) 

(A3c) 

(A4a) 
(A4b) 
(A4c) 
(A4d) 



APPENDIX B: THREE-PHOTON ASYMPTOTIC OUTPUT STATE 

In this Appendix, we present the asymptotic output state after scattering a three-photon right-going Fock state off a 3LS or 
4LS. The form of the wave functions is 

f dkidk2dk^-^a(kx)a(k2)a{ki)\cp^^Hkx,k2M)), 
J V3! 

dxldx2dx3[—tttk^^l,^^k^{xl,X2,X3)a\{xl)al{x2)al{xi) + —ttrk^ 

+ j^trrkt,k2M(^i'~^2,-X3)a]^{xi)al{x2)al{xi) + ^rrric^ . (Bl) 

Here, ttti,j,k2.k}ixi7X2,x-i), ff?"i:(_<;2,i:3(jci,X2,X3), trrkj,ki.k}{xi,X2,X3), and rrrin^k2,k3(xi,X2,X3) are the terms representing three- 
photons being transmitted, two being transmitted and one reflected, one being transmitted and two reflected, and all three being 
reflected, respectively. They take the following general form {a,p,y -t orr) 

a/3yk^,k2,k3(^i,X2,X3) = Y^akQ^/3kQ^ykQ^hkQ^{xi)hkQ^{x2)hkQ^{x3) + jY^[akQjikQ^(xi)Bf^ , (x2,X3) 



+likQ^ hkQ^ i^2)Bfl. (XI ,X3) + ykQ^ hkQ^ (^3)Bg (^1 ,^2)] + \Yj^\ 

PQ 



(3) 

P, ."^2 '"''3 ' 



(B2) 



where and are the single-photon transmission and reflection probabilities given in Eq. (fT3b). ^[^'^^ (■'^i ' -^2) is given in Eq. ( fTSl ). 
and i^^{x\,X2,X3) is given in Eq. ( fTOb . In Eq. ( IB2l i. the first term comes from the process of three-photons passing by the 
atom as independent particles. The second term corresponds to the process of one-photon passing through as an independent 
particle while the other two photons form a composite particle in a two-photon bound-state (with three possible combinations). 
The third term originates from the three-photon bound-state process. 

APPENDIX C: SECOND-ORDER CORRELATION FUNCTION IN THE SCHRODINGER PICTURE 



In this Appendix, we demonstrate the equivalence between Eq 
correlation function is defined in the Heisenberg picture as, 

{4f()\a\xi, t\)a\x2, t2)a(x2, t2)a(xi , f 1 )|i/ro) 



and Eq. ( l27l i in the main text. Typically, the second-order 



g^^\xi,ti;x2,t2) = 



{ifro\aHxi,ti)a{xi,ti)\4ro}{iffo\a^X2,t2)a{x2,t2Mo} 



(CI) 
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where |i^o) is the state in the Heisenberg picture, or equivalently, the initial state in the Schrodinger picture and a^{x,t) is the 
operator in the Heisenberg picture, a^x, t) can be expressed in terms of the operator in the Schrodinger picture as 

a\x,t)^e'"'l^' aHx)e-'"'l\ (C2) 

Taking xi - X2 - x m Eq. dClb . we obtain the two-time correlation function 

(2), . . s {tlJo\a^{x,ti)a^{x,t2)a{x,t2)a{x,ti)\ipo) 

'{x,ti;x,t2) 1 ; ■ (C3) 

{tffQ\aT(x, ti)a{x, f i)|i/'o)<i/ro|aT(x, t2)a{x, t2Mo} 

If the field operator satisfies the following relation 

a\x,t)^aHx-ct), (C4) 

g^^\x,ti;x,t2) is then the same as g'-^\x,ti;x' Ji) with x' - x — c{t2 -ti). Using Eqs. iC2\ and (IC4b . we can rewrite ( IC3I ) in the 
Schrodinger picture as 

(2) f mti)\aHx)a\x')a{x')a{xmti)} 

where \tl/(ti)) is the state at f = fi evolving from the initial state li/^o) under the Hamiltonian H. Therefore, as long as Eq. ( |C4| ) 
holds, the definition of g^^^ in the Heisenberg picture Eq. iC3i is equivalent to Eq. iC5i defined in the Schrodinger picture. 
Physically, this means that measuring the two-time correlation at the same spatial position is equivalent to measuring the spatial 
correlation at the same time for a non-dispersive field. 

In our problem, it is straightforward to show that Eq. iC4i is satisfied by the right-going field. With the Hamiltonian defined 
in Eq. (HJ in the main text, the equation of motion for the right-going field in the 4LS case is 

(yx^-cj!) " T ^"2(0 + S UD] 6ix). (C6) 

Formally, the above equation can be integrated to yield 

alix, t) = a\j,.Jx - cf) + ^ [S - xjc) + S " ^Ic)] 6{x). (C7) 

A similar expression can be obtained in the 3LS case. Hence, Eq. ( IC4b holds, and we use Eq. ( IC5b to evaluate the second-order 
correlation function of the transmitted field with |i/'(fi)) being our final output state. 

I 
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